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In the standard geometric approach, the entanglement of a pure state is sin 9, where 6 is the 
angle between the entangled state and the closest separable state of products of normalised qubit 
states. We consider here a generalisation of this notion by considering separable states that consist 
of products of unnormalised states of different dimension. The distance between the target entangled 
state and the closest unnormalised product state can be interpreted as a measure of the entanglement 
of the target state. The components of the closest product state and its norm have an interpretation 
. in terms of, respectively, the eigenvectors and eigenvalues of the reduced density matrices arising in 

^ l' the Schmidt decomposition of the state vector. For several cases where the target state has a large 

degree of symmetry, we solve the system of equations analytically, and look specifically at the limit 
(~| , where the number of qubits is large. 

^ ' PACS numbers: 03.65.Ud, 03.67.Mn 
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■ With recognition of its role as a resource in quantum computing the nature of entanglement in quantum systems 
' is a problem of much current interest Of particular importance is a quantitative measure of entanglement j5|. 
, Two of the more commonly used measures are the von Neumann entropy, which is based on reduced density matrices 

and a geometric measure, which is based on the distance to the nearest product state [6'-^10|. In this paper we 
, introduce a geometric measure of entanglement based on the distance between an unnormalised product state and a 

■ target entangled state. The norm of the closest product state can be related to both the distance and angle between 
^ I the product and target states. This result motivates the interpretation of the distance to the closest product state as 

a measure of the entanglement of the initial state. 

We begin by defining our notation. We consider a system of q qubits. The dimension of the corresponding Hilbert 
space is n = 2"^. We will decompose the system into a set of subsystems. The subsystems are labelled A, B,C, ■ ■ ■ 
They have dimension u,v,w, - ■ ■ such that n — u ■ v ■ w ■ ■ ■. An arbitrary set of basis states of system A is labelled 
the basis states of B are the basis states of C are \k), etc. Using this notation we write: 

It— 1 v — l w — 1 

\A) = Y,a,\^), \B)^Y.b,\j), |C) = ^ c,:|fc), ... (1) 

i=0 j=0 k=0 

We consider an arbitrary normalised entangled pure state and write its wave-function: 

1^) = E E E • • • ^'^■'=•••1*) ® 1^') ® • ■ ■ ' (v^iv^) = 1 ■ (2) 

i^O j^O k^Q 



- 
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In this paper we introduce a new geometric measure of the entanglement of this state. The paper is organised as 
follows. In section In] we introduce our geometric measure of entanglement. In section Hill we show that, for a given 
entangled state, a connection can be established between the components and norm of the closest product state, and 
the basis states and eigenvalues of the Schmidt decomposition of the entangled state. In section IIVI we study some 
general symmetries of our measure. In section IVl we derive some exact solutions for cases where the target state 
has a large degree of symmetry, and in section IVTl we present our conclusions. 



II. OPTIMUM EUCLIDEAN DISTANCE 



In this section we introduce a new geometric measure of entanglement. We look at the distance between the pure 
entangled state ^ and an arbitrary unnormalised product state. Extremizing this distance allows us to identify the 
closest product state. The distance between the state and this closest product state is our geometric measure of 
the entanglement of 

We consider the product state: 

u—l v — 1 w — 1 

\(t>) = \A) (g) |S) (g) |C) «) . . . = ^a,\i) <^^bj\j) «) ^ Cj\k) (g) . . . (3) 

i=0 j=0 k=0 

The state \(j)) is not assumed to be normalised: 



v-l 



{c^\c^} ^NaNbNc...; iVA = (A|A) =^a*a„ A^s = - ^ ••• (4) 

1=0 j=0 



The distance between the states and is: 



= 1 - (01^) - (tAI^) + TVaA^bA^c • ■ • 

' u— 1 V — 1 W — 1 \ 

XI XI XI ■ ■ ■ • ■ • - X*jk-) {atbjCk ... - Xijk-) ■ (5) 



1=0 j=Q k=0 

We extremize this distance with respect to the coordinates of \B 



dai 
dck 



- ^ a*NBNc • ■ • = E E • • ■ I ^J^fe • • -^-fc.. 

/u—l w — 1 

= 6*iV^7Vc...= EE- 



O-iCk ■ ■ ■ Xijk... 



yi=0 k=0 
u—l v—1 



= cINaNb...= I EE"- I '^rb,...x*,k... 

i=0 j=0 



Rewriting Eq. ©, we have: 



' u—l v—1 w—1 



(6) 



NaNbNc •■•- EEE"- • • • ^^k... = im ■ (7) 

^4=0 j=0 k=0 J 

In exactly the same way we could extremize the distance in Eq. ([5|) with respect to the coordinates of {A\, {B\, . . . and 
obtain: NaNbNc . . . — ((/'IV')- Substituting these results into ([5]) we find that at the extrema the distance between 
It/;) and \(f>) is: 

= 1 - NaNbNc . . . = 1 - cos^ 6c , (8) 
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•J{<Pn\'I>n) = 1 



FIG. 1: Geometrical comparison of the entanglement measure using unnormalised {(f}) and normalised (<^]v) separable states. 



where we have defined the critical angle 0c as the angle between \tp) and |(/)) at the extrema: 



cos Uc 



^NaNbNc . . 



(9) 



critical 



In order to demonstrate the consistency of these results, we look at the Cauchy-Schwartz inequality which requires: 



Using (jll and ([5]) we have: 



. {NaNbNc ■■■)> \{NaNbNc ■■■)- {'^\4>)Y 



(10) 

(11) 

(12) 



Using ([7]) and (jS]) we find that at the extrema 

1 - NaNbNc • • • > ^ NaNbNc • • • < 1 , 

This inequality guarantees that the cosine of the critical angle (Eq. ([9])) is a real number which satisfies 
1 > cos^c > — 1: and that the square of the critical distance (Eq. ([5])) is real and positive. 

We can compare the results in Eqs. ([8]), (jl]) with those that would be obtained using a product of normalised states. 
We use = {B\B) =■■■ = !. We take the derivative of with respect to as before, but now we insert a 

Lagrange multiplier term of the form A(A|A). This approach produces the same critical angle as before (Eq. ©). 
The corresponding minimal distance is: 

2 



D 



N 



/'iv - V'l'^JV - ^) = + (l - v/(M) -2(l-cos0c): 



(13) 



which can be compared with the result for unnormalised product states (Eq. ([8])). 



III. CONNECTION TO THE SCHMIDT DECOMPOSITION 



The equations in ^ which determine the extremal points of the distance to the closest product state are non-linear 
and must be solved numerically, except in special cases. One of the special cases for which a closed-form solution exists 
occurs when the n-dimensional system is decomposed into two subsystems. We consider a u-dimensional subsystem 
A, and a w-dimensional subsystem _B, such that n = uv. In this case the equations ^ decouple to yield: 

ti — 1 v~l 

aiNANB = ^^XtjXi'jat' , 

v — 1 u — 1 

bjNANB = EE^^^x:,'^/- (14) 

j'=0 i=0 



Each of these equations can be solved for the product NaNb and give, respectively: 

v-l 



det 



det 



u-l 



= 0. 



0. 



(15) 



These solutions can be related to the eigenvalues of the reduced density matrix. We consider the n-dimensional 
state in Eq. ^ and write its wave-function and density matrix in the computational basis: 



n— 1 n — 1 



(16) 



Z = 2 = 2'=0 

Decomposing the system into a u-dimensional subsystem A and a u-dimensional subsystem B. we obtain: 

1^) = E E x^^i*) ® 1^') ' - E E E E xnxir\^^'\ ® \3){j'\ ■ 

1=0 j=0 i=0 i'=0 j=0 j'=0 



(17) 



Next we calculate the reduced density matrices. The reduced density matrix pA is obtained by tracing over the 
subsystem B, and the reduced density matrix ps is obtained by tracing over the subsystem A. The definitions are: 



u-l 



TlAip) = ^((i|®li3)p(|t)®lB) , 
t=0 

■u-l 

Trs(p) - ^(1a® (t|)p(lA®K)) , 



(18) 



t=0 

where 1^ and 1 b are the identity matrices in the subspaces of A and i?, respectively. We obtain: 

v—l u~l u—1 v—1 v — 1 v — 1 u — 1 u—1 

PA = Trsip) ^J2Y.T.J2Y. x^,xh'm{f\m{^'\ = E E E x..x:',n)(^'i , 

t=0 i=0 i'=0 j=0 j'=0 j=0 i=0 i'=0 



(19) 



or, in terms of components: 



v-l 

iPA)^^' = E ^vXl'j 
3=0 



(20) 



u-l 
i=0 



Similarly, the reduced density matrix ps = Tryi(p) is obtained by tracing over the subsystem A and can be written: 

(21) 

Using Eqs. ((20|) and ([2T|) we can rewrite the extremal conditions of Eqs. (|T4l) in the form: 

(/3A)ii/ ai> = NANsai , 

A^^iYsfe, . (22) 



i=0 

D-1 



E ^PB)n' 

Equation shows that NaNb are the eigenvalues corresponding to the eigenvectors Ui and 6j of the reduced density 
matrices pA and pB- 
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This result can be interpreted in terms of the Schmidt decomposition as follows P, [ll| - [lH . Let us write the extrema 
conditions of Eq. ^ as 



M-l 

bja ^ ^a*Xzj (23) 



where — ai/ \/ Na and bj — bj / \J Nb are, respectively, the left-singular and right-singular vectors corresponding 
to the singular values a — ^NaNb of the matrix x- The singular value decomposition of the matrix x can then be 
written as 



X 



AJ:B^ (24) 



where the columns of the (unitary) matrices A and B are, respectively, the vectors ai and bj, and I] is a diagonal 
matrix whose elements are the singular values ct. This can be used to rewrite the state j?/;) in ^T7\ in terms of the 
Schmidt decomposition involving a single summation: 

min(u— l.ij — 1) 

IV') = E Vn\ak)'E>m, (25) 

fe=0 

where the Schmidt coefficients pk, which are identified with the singular values a, satisfy X^fcPfc ~ 1' ^^^1 the states 
\ak) and |/3fc) are identified with, respectively, the left-singular and right-singular vectors and bj. Calculating the 
corresponding reduced density matrices using (|18p we obtain: 



PA = TrB{p) ^^Pk \ak){ak\, 

k 

PB = Tr^(p)-^p, \f3k}{f3k\, (26) 

k 

from which one can see that \ak) are the eigenvectors of pA ~ Ttb{p) and \/3k) are the eigenvectors of pB — Tr a{p) 
with corresponding eigenvalues pk- 

We remark that the Schmidt decomposition gives rise to several other measurements of entanglement. The reduced 
density matrices pA and pB in p6p have the same non-zero eigenvalues and, for a product state, only one non-zero 
eigenvalue is present. As a result, for a non-product state, the eigenvalues of the reduced density matrix can be used 
to quantify the degree of entanglement. One commonly used measure of entanglement is the von Neumann entropy: 

S = -Tr (pa In pa) = -Tr (ps Inps) = - ^Pfc log2Pfe • (27) 

k 

It is interesting to consider the particular case that the n-dimensional space is split into the product of a single 
qubit space and another space of dimension u — njl. In this case, one of the equations in (|15l) will become a quadratic 
equation for the product NaNb with solutions: 



NaNb = M± = 2 



l±yT^4Cj, C = J2J2\^OjXik-XijXokf ■ (28) 

j=l k=0 



The Schmidt decomposition (Eq. p5t ) becomes: 

1^) = V7ZT|a+) ® 1/3+) + ^/7^ |a-) ® I/3-) . (29) 

Using Eq. ([9]) we can relate the eigenvalues NaNb — p± to the cosine of the angle between and the closest product 
state \ip). Since /i+ is the larger of the two eigenvalues we write cos6'„iax = p+ and Eq. becomes: 

IV-) = cos6l,„ax |a+) (E) 1/3+) + sin6'n,ax \a-) (E) |/3_) . (30) 

If 1-0) was a product state, the reduced density matrix would have only one non-zero eigenvalue and we would have 
/i_ = sin'^ (?max = 0. This result is consistent with our geometric approach: from Eq. ([5]), sin^ ^max = corresponds 
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to a zero minimal distance between the state \ip) and the nearest product state, which means that lijj) is itself a 
product state. 

The Schmidt decomposition can also be applied to multipartite pure states [T^ . One starts with a state \ip) and 
decomposes it into two subsystems: a single qubit system A, and a subsystem (BC . . . Z) containing all other qubits. 
Using a Schmidt decomposition we can write: 



i,BC...Z\ 



(31) 



One then decomposes into two subsystems: another single qubit system B, and a subsystem {CD ...Z) 

containing all other qubits. Using a Schmidt decomposition we can write: 



i^r-^) - E « m^-') • (32) 



This process is continued until the last two qubit spaces Y and Z are reached, with the result: 



E 



Now we consider the geometric interpretation of this result. Consider the distance between \tp) and a state 
|0A;_BC...z^ _ i^BC.z^^ where 10"^) is a single qubit state and l^^*-^ - -^) is the state of the remaining qubits. 
This distance is: 

D%Bc...z = (V^ - [0^ ® ^^^-1 k - [<t>^ ® 0^^-^] ) . (34) 

Finding the extremal points of this distance will result in a system of (linear) equations, as in Eq. (jl5p . which 
determine the components of the state - There is a direct correspondence between the coefficients of the 

closest product state and the basis states of the Schmidt decomposition (Eq. (pijl ). and the norm of the closest 
product state is related to the corresponding Schmidt coefficients. 

We can then consider the distance between the state and a state |0S;CD...z^ ^ j^s^ ^ ^^CD...zy^ 

D%CD...z = ('^^^■■■^ - [0^ ® <^^^-1 k^*"-^ - [^^ ® 'f''-^\ ) , (35) 

where is a single qubit state and \(^'-^^---^) is the state of the remaining qubits. Extremizing this distance 

will again result in a system of linear equations determining the components of the state |(/)^''^^ - '^). Once again, 
there is a direct correspondence between the coefficients of the closest product state and the basis states of the 
Schmidt decomposition (Eq. ([5^ ). and the norm of the closest product state is related to the corresponding Schmidt 
coefficients. 

This process may be continued until the last two qubit states |(/)^) and are reached. In analogy with Eq. ([3]), 
we can define the cosine of the critical angle Qc'- 



COS Oc = 



... 10) = (36) 

critical 



The quantity sin^ 6c = I — cos^ 9c is a measure of the entanglement of the original multipartite state. 

The advantage of this procedure is that it involves solving a series of linear equations, as compared to the approach 
of Section [IT] which produces the non-linear equations of Eq. ([6]) . The disadvantage is that the result depends on the 
order that the series of decompositions is made: the sequence ABC . . . YZ described above will difi'er from the sequence 
BC...ZA. In Ref. the entanglement is given by the minimal value obtained by looking at all permutations of 
the possible orders of the decompositions. 



IV. SYMMETRIES 



The equation that gives the distance between the target entangled state and the nearest product state (Eq. ([8|)) is 
invariant under certain transformations of the parameters {a^, bi, . . .}. In order to see these symmetries explicitly we 
can reparameterize each set of coefficients. For the coefficients we write: 

a, = A,e*<. (37) 
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with similar equations for the coefficients {bi,Ci, . . .}. Using generahsed spherical coordinates in n dimensions we can 
rewrite the set of real variables Ai in terms of a magnitude A and (n — 1) angles 0;: 







A2 




A3 




Ak 


= A 


An-l 




\ A„ J 





I 



COSf^i 

sin Qx COS 6*2 
sin Q\ sin Qi cos Qi, 



\ 



(nti'sine.) 



cos Vk 



smf^i • • • sm6'„_2 cosyn-i 



We can also parameterise the phase angles a" as: 









1^ 1 ^ 




/ 


g^k ^ 








gi(Q!2-ai) 
















EE e*"i 












gi(Q„-Ql) ^ 




v 


g.;3„ ^ 



where /3i — 0. Using this notation we write: 

We can use this parameterisation to rewrite the distance function ([8]). We make the definitions: 

^N\^), (010) = 1. 

Note that the first line in (HTj) gives: 

{^\cb) ^A^B^C^-^-^NANsN^- 
'Usmg these definitions the distance in Eq. ([SJ becomes: 



= -N 



im + (01^) 



1. 



Extremizing we obtain: 



~dN 



^ 2N 



(^10) + {4>\ij) 



TV = Re 



which means that at the extrema: 



(V'I'^)J , 

Dl^N^ - N{2N) + 1 = 1-N^ . 



(38) 



(39) 



(40) 

(41) 
(42) 

(43) 



(44) 
(45) 



We note that Eqs. dH]) and (US]) are consistent with 

To make this result more clear, we can look explicitly at the dependence of the distance function on the overall 
phase of the coefficients of the product state. We define the overall phase angle S — + a\ + . . . and obtain: 



(01^-) = Ni^lij) = Ne-^' hine-'"^ fj{0')e-'^'^ ■ ■ ■ X^r 

i=i j=i 

U V 

im = ^^(^-10) = Ne^' EE"- Mne'^^ fji&'y^'^ ■ ■ ■ x*,.- ^ 

i=i j=i 



(46) 
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which gives: 



— (<^|^) = N—{^\^p)^-tN{m. (47) 



(V-l^/-) + (01V) =-«^ (V|<^) - (<^|V) =2mm{m, (48) 



From Eq. we have: 

dS dS 

which means that at the extrema, lm{'ip\(j)) — 0. From Eq. (1441) we obtain that at the extrema: 

N = (VI0) = (01V) . (49) 

This means that the critical angle defined by: 



COS Oc = 



is the same critical angle as in Eq. ([9]) 



= N, (50) 

critical 



V. EXACT SOLUTIONS 



In this section we look at some states with a large degree of symmetry for which the equations ([6]) can be solved 
exactly. We consider a system of q qubits and divide the Hilbert space of dimension n = 2^ into q spaces, each of 
dimension 2. Using the notation of section U we have n = u ■ v ■ w ■ ■ ■ with u = v = w — ■■■ = 2. The basis states in 
each single qubit system are: 



1\ /O 




, |1) = 1 ■ (51) 



The product state in Eq. ^ becomes: 



(ao|0) + ai|l))(6o|0) + 6i|l»(co|0) + ci|l)) • • • 

(aoboco ■ ■ ■)\0, 0, 0, • • •) + (ao6oCi • • OlO, 0, 1, • • •) + {aohci • • OlO, 1, 1, • • •) + • • • (52) 



We look for a solution of the form: 



which means Na = Nb = ■ ■ ■ = N. 



do = bo = ■ ■ ■ = ao , 

ai = 6i = • • • = ai , (53) 



Case 1 

Consider: 

IV) = VpIOO ...0) + yr^|ll ...1). (54) 
The only two non-zero components of ipijk--- are: 

Voo---o = Vp, 

Vii.--i = a/1 -p, (55) 



and Eq. © becomes: 



Using Eq. Eq. becomes: 



aaNsNc ■ ■ ■ = y/p ^oCo • • • 
aiNsNc--- = v^l -pbici 



Solving this set of equations we obtain: 



al 



which gives: 



Rearranging we obtain: 



N = al + a\ = 



J^2(q~l)/(q~2) 
(1 -p)l/(9-2) ' 

pl/(g-2) ' 

7V2(g-l)/(g-2) 
(1 -p)l/(<?-2) 

p{l-p) 



(l_p)V(^ 
^ pl/(9-2) 



[pl/(9-2) + (1 _p)l/(«-2)j9 2- 

We can look at the large q limit. Defining A = [^1/(9^2) ^ _ ^-j 1/(9-2) j "J ^ jja,ve, 
log A = ((j-2)logU/(«-2) + (f-p)i/(9-2)l ^ 



(g-2)log(pi/(^-2))^(^_2)log 



log p + (g - 2) log 2 + i log ( 

2 \ p 

(g-2)log2+ilog[p(l-p)] , 



l~p 



1/(9-2)1 



log 2'-Vp(1-p) 



which gives 



TV? 



v/p(l-p) 

29-2 



Consider: 



Case 2 



|^) = [|fO...O)+ |0f ...0> + ... + |00 ... f)]/V9. 



The only non-zero coefficients ipijk--- are: 



'01OO---O — V-'oio---o — V'ooi---o — • • • — V'ooo-'-i — —;= , 

V9 



and Eq. © becomes: 
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qoNbNc ■■ ■ = [bicodo h 6oCido H feoco^i H • ■ •] 

aiNsNc ■ ■ ■ = — fcoco^o ■ ■ • (65) 



Using Eq. ((53| . Eq. ((65ll becomes: 



ai7V9-i = ipal~\ (66) 



Solving the set of equations in (I66p we obtain: 



1 (g- l)(9-i)/2 ' 



which gives: 



^2(9-l)/(<;-2) 

(g-l)/(g-2) 



^ = "o + «N— t-tttt;-^- (68) 



1- i 
9 



Rearranging we obtain: 



IN 9-1 



7V« = ( 1 - - ) . (69) 



Case 3 



Consider: 

IV') = [|1 100 ...0) + 101 10 ...0) + 1001 1 ...0) + ... + |00 ... 11) + |10 ... 01)]/vg. (70) 
The only non-zero coefficients ipijk--- are: 

V'llOOO---0 = "^OllOO-'-O — ■0OO11O---O — ^00011---0 — ■ ■ ■ = ■0OOOOO---11 — ■01OOOO---O1 = —pi ■ (71) 
Using Eq. dSS]), Eq. (O gives: 

iV-i = ^alar' - ^alar' . (72) 
2' 

^r(9-i)/(<z-2) ^ (73) 



Solving the equations in d72|) we obtain: 

ao 



Qi -'2 
1 1 



which gives: 



22-9 „2{q~2)' 



N ^al + al = ^_2i/(2(2-9))^i+i/(g-2) ^2(g-i)/(,-2) _ (74) 
2 

Rearranging we extract: 

7V« = 4(g - 2)«^2^i-« . (75) 
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Case 4 

Consider the case where has q qubits consisting of all possible combinations of p entries of "1" and q — p entries 
of "0". The normalisation is: 



r 



(76) 



Using Eq. Eq. ^ becomes: 



9-1 



(77) 



Solving the set of equations in ([77|) we obtain: 



,,9-2 _ 



AT?- 



an = «! 



2 / 9 



(78) 



which gives: 



N = al + al= ( - 



(p-l)/(g-2) 



^2(q-l)/(g-2) 



1 - £ 
9 



(9-p)/(9-2) 



1/(9-2) 



(79) 



Rearranging we obtain: 



1-P 



P~i 



(80) 



In the limit q ^ p limit, we can approximate: 



q-l\_ (g-1)! 



p-ij ip~iy.{q~py. {p-iy. 

« exp(-p) , 



^—qP-^+0{q^-') 



q 



(81) 



which leads to: 



p'^^cxp(-p) /I 



(82) 



A. Results 



We show below a graph of our results for the entanglement measure as a function of the number of qubits for the 
four cases solved in this section. 
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FIG. 2: The entanglement as a function of q from Eq. ([8)l and Eqs. (|62p - dotted/red, (|69l) - solid/green, ((75} - dashed/blue 
and H82p - dot-dashed/magenta. 



VI. CONCLUSIONS 



We have considered a generalisation of the usual geometric measure of entanglement of pure states using the distance 
to the nearest unnormalised product state. This definition does not lead to any computational advantages, since the 
set of equations that determine the measure are still non-linear in general. However, our definition does provide an 
interpretation of the standard entanglement measure as the distance to the closest product state. We have also found 
a relationship between the norm and components of the closest separable state, and the coefficients and basis states of 
the Schmidt decomposition of the state For several cases where the target state has a large degree of symmetry, 
we have solved the system of non-linear equations analytically, and looked specifically at the limit where the number 
of qubits is large. These results indicate that our new definition of entanglement, while similar to other definitions 
that can be found in the literature, is worthy of further study. 
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